Inducing odd-frequency triplet superconducting correlations in a normal metal 
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This work discusses theoretically the interplay between the superconducting and ferromagnetic 
proximity effects, in a diffusive normal metal strip in contact with a superconductor and a non- 
uniformly magnetized ferromagnetic insulator. The quasiparticle density of states of the normal 
metal shows clear qualitative signatures of triplet correlations with spin one (TCSl). When one 
goes away from the superconducting contact, TCSl focus at zero energy under the form of a peak 
surrounded by dips, which show a typical spatial scaling behavior. This effect can coexist with a 
focusing of singlet correlations and triplet correlations with spin zero at finite but subgap energies. 
The simultaneous observation of both effects would enable an unambigous characterization of TCSl. 

PACS numbers: 73.23.-b, 74.45. +c, 85.75.-d 



Hybrid superconducting/ferromagnetic circuits allow 
to obtain fascinating unconventional superconducting 
correlations [Ij. Although a standard s-wave supercon- 
ductor naturally hosts even-frequency superconducting 
correlations, the superconducting proximity effect in an 
adjacent ferromagnet can lead to odd-frequency pairing, 
because the ferromagnetic exchange field lifts time re- 
versal symmetry. This can correspond to either triplet 
superconducting correlations with spin zero (TCSO), i.e. 
correlations between opposite electronic spins, or triplet 
superconducting correlations with spin one (TCSl), i.e. 
correlations between equal spins. The existence of TCSO 
is confirmed experimentally since a decade (see e.g. 0- 
0]). The possibility of obtaining TCSl has been inves- 
tigated more recently [H']. The features observed so far 
are quantitative. It has been observed that, in certain 
conditions, ferromagnets can sustain a supercurrent on a 
much longer lengthscale than expected 'i']. This suggests 
the presence of TCSl, because in diffusive ferromagnets, 
TCSl are expected to propagate on a much longer dis- 
tance than TCSOpT^. Most of the strategies discussed so 
far to observe TCSl require to measure a supercurrent, 
which is an energy-integrated quantity. Alternatively, the 
quasiparticle density of states (DOS) yields spectroscopic 
information and thus appears as a powerful tool to char- 
acterize the nature of superconducting correlations. This 
paper presents a geometry in which the DOS gives a par- 
ticularly rich access to odd-frequency superconducting 
correlations. 

I suggest to use a lateral geometry, where the supercon- 
ducting correlations are induced in a normal metal (NM) 
strip in contact with a superconductor (S) and a ferro- 
magnetic insulator (Fl) with two non-colinear magneti- 
zation domains. This configuration is compatible with 
spatially resolved DOS measurements using several tun- 
nel contact probes or a low temperature STM[^. The 
propagation of odd-frequency superconducting correla- 
tions has been studied thoroughly in ferromagnets [![, but 




FIG. 1: Scheme of the lateral geometry considered in this 
work (see text). 



only elusively in NMs lol- 1^ . Here, the Fl turns the 
NM strip into an effective ferromagnet with an unusu- 
ally weak exchange field. This is advantageous to study 
the propagation of odd-frequency superconducting corre- 
lations. The TCSl induce a zero-energy peak in the DOS 
of the NM0. Such an effect is not specific to TCSl[ll- 
16|. However, in the present geometry, important ad- 
ditional features allow to identity unambiguously TCSl. 
Indeed, the low-energy DOS peak is surrounded by dips, 
and this structure shows a characteristic spatial scaling 
behavior when one goes away from the superconducting 
contact, i.e. it "shrinks" on a lenghtscale which depends 
on a longitudinal Thouless energy. If one observes simul- 
taneously finite energy dips, which confirm the existence 
of the effective exchange field inside the NM, the zero- 
energy peak points unambiguously to TCSl. 

I consider the lateral geometry of Fig.l, where the cen- 
tral element is a NM strip with thickness /i, width w, and 
longitudinal coordinate x. A portion with length (11 of 
the NM strip is contacted to a S and a Fl domain magne- 
tized along a direction rfiL for a; < 0, and a portion with 
length dji is contacted to a Fl domain magnetized along 
a direction for a; > 0. I assume that the structure is 
diffusive, so that a quasiclassical isotropic Green's func- 
tion G can be used to describe the propagation of super- 
conducting correlations inside the NMjl9|. The function 
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Y-^0.005, Ej =0.5A, E^=10A, r=0.05A 
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FIG. 2: Absolute values of the anomalous components fs, ft 
and ft of G in the NM strip at e = (left panel) and e = Et^^' 
(right panel) respectively, for m_B = z and dj^j^) — >■ +oo. 



G has a structure in the spin and Nambu (electron-hole) 
subspaces, so that it can be decomposed in terms of the 
spin[Nambu] Pauli matrices ai[Ti] (see below). In the fol- 
lowing, I use rfiL = X. In the colinear case mn — x, only 
singlet correlations and TCSO can appear inside the NM. 
In the non-colinear case m^j ^ x, TCSl can also appear. 

The spatial evolution of G inside the NM is described 
by the Usadel equation hPpy.j — —{is — T) [f^, G], with 
p the resistivity of the NM[20| . The density of matrix cur- 
rent j = GSIG/ p characterizes the flows of charge, spin 
and electron-hole coherence in the device. The rate T ac- 
counts for inelastic processes. The Usadel equation alone 
is not sufficient to predict the behavior of G, because the 
influence of the S and FI contacts must be taken into 
account. When h and w are smaller than the typical 
spatial scale ^ characterizing the variations of G inside 
the NM 2l[ , it is possible to derive a one-dimensional ef- 
fective Usadel equation on G ~ G{x), which takes into 
account the effects of the S and FI contacts. This equa- 



tion can be derived with an approach inspired from cir- 
cuit theory [2^, using spin dependent boundary condi- 
tions for isotropic Green's functions [23*]. This requires to 
define a surface tunnel conductance Gfp for the S/NM in- 
terface and a surface conductance-like coefficient G^ for 
the NM/FI interface. The last parameter accounts for 
the fact that electrons from the NM are reflected by the 
FI with spin-dependent reflection phases. Indeed, the 
internal Stoner exchange field of the FI can affect the 
evanescent tails of electronic wavefunctions on a scale of 
a few atomic layers. One finds, for a; < 

2/2V,(GV,G) = -27,[f3,G]-7T[G,G5]+i70[T3fTL,G] 

(1) 

with ETh = hD/h'', 7e = {ie - T)/EThnT = G?rph, 
70 — G''^ph and o'i(fl) — mLj^jiyd. Here, Gs denotes the 
value of the isotropic Green's function inside the super- 
conductor. For X > 0, one finds a similar equation with 
7t replaced by and by crj^. 

Equation ([1]) describes an interplay between the super- 
conducting and ferromagnetic proximity effects. On the 
one hand, the 7t terms tend to induce a minigap inside 
the NM, due to the confinement of electrons in the z di- 
rection. This effect depends on the interface parameter 



7t but also the lateral Thouless energy Eth defined be- 
low Eq. ([!]). For instance, at the left end of the NM strip, 
the DOS is suppressed for energies |£| < A = EThlr/^ 
when ^ -|-oo, EthIt A and 7^ = 0. On the other 
hand, the interface parameter G^ causes an effective ex- 
change field E^JJ = EThl(i>l'^ — fiDG^p/h oriented along 
X [rhu] at the left[right] side of the NM strip, i.e. a; < 
[x > 0] [i^l- This type of effective field has already been 
observed and accurately characterized for various types 
of S/FI bilayersjii-ii]. For a; < 0, the exchange field 
due to the left contact splits the minigap induced by the 
S contact and induces TCSO with respect to the x di- 
rection. These correlations can propagate to the right 
part of the NM strip where they correspond to TCSl 
when mR = z. The induction of TCSl in diffusive S/F 
structures present similarities with this process [Hj. How- 
ever, minigap effects are usually destroyed in ferromag- 
nets, due to strong ferromagnetic exchange fields, except 
in some particular disordered cases (2^. Besides, E^l^ 
is expected to be much smaller than the exchange field 
inside standard ferromagnets. Hence, the situation stud- 
ied here is qualitatively different from that of standard 
S/F structures. Interestingly, for a given type of NM/FI 
contact, the amplitude of E^l^ can be tuned by choosing 
h at the sarnple fabrication stage, since E^^-f scales with 
26l 30l | . This gives an interesting fiexibility with 
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respect to material constraints. 

In the following, I assume that the value G5 of G inside 
the superconductor is equal to the bulk BCS value, i.e. 
Gs = cos{9s)f3 + sin(0s)fi, with 63 — arctan[A/(— ie -|~ 
r)]. Inside the NM, one can use the angular parametriza- 
tion G = T3(coshry) cos(6') +iamh{i]) sin{9)v.(T)+Ti{fs + 
ft.ff), with fs = cosh(?7)sin(6'), ft = -i sinh(?7) cos(6')t/, 
and V a unit vector. This convention automatically sat- 
isfies the normalization condition G^ = 1. I discuss be- 
low the anomalous components of G, i.e. fs, ft — ft-x 
and fl = ft-z, which reveal the existence of supercon- 
ducting correlations inside the NM. Defining TCSO and 
TCSl components requires to define a reference direction 
rhref- A natural choice is to use rhref = tol for a; < 
and rriref = ttt-r for a: > 0. Then, TCSO and TCSl corre- 
spond to the components of ft parallel and perpendicular 
to rhref, respectively. The DOS inside the strip can be 
calculated as N{e)/No = Y.a R-e[cos(6'-Kicr77)]/2, with iVo 
the normal state DOS of the NM. 

The spatial evolution of the DOS can be first studied 
with an analytic linearized description, which yields a 
transparent interpretation of the circuit behavior. This 
approach is valid in the limit of a weak superconducting 
proximity effect, i.e. 0, 77 ^ 1 inside the NM. This occurs 
e.g. when 77- is small and F large, so that the minigap is 
closed and just gives residual dips at the left side of the 
NM strip. For simplicity, I assume dj^i^n^ — > -\-oo. For 
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FIG. 3: Density of states in the NM in the non-colinear case 
ruR = z (panel a), the colinear case ifiR — x (panel b), and 
the case where there is no FI contact at a; > (panel c) . The 
data are shown in terms of the deviation 8N = N — No of 
the DOS in the NM from the normal state value A^o- The 
parameters used are the same as in Fig. 2. 



< one has 



[/„ , = [0, 0, l]Aoe'=-^+ ^ [1, a, 0] + ^.e 



a=±l 



e 



with -7tCOs(6's) = -2i[{e + iT)/ ETh], and fc^/i^ 
7T 008(6*5) = 2i[<TEllf - £ - iT]/ETh- Here, I note 



(2) 



8^ = arctan[ 



7t sin [0s] 



7T cos [0s] + 2[-ie + r + iaE'JJ]/ETh 



(3) 

the value of — ii^a at a; — oo. For x > Q and rfiR = 
one has 



[/.,/f,/t1 = [0,l,0]i?, 



be 



— A,'o X 



[l,0,a]B„ 



(4) 



with k'^h? 



= -2i[{e + iT)/ETh] and k^h^ = '2i[(TEgf ~ 
e — iT]/ETh- For x > and rnji — x, the second and 
third components of the vectors in Eq.(j4]) must be ex- 
changed. The coefficients A^, Aq, Ba, and Bq can be 
calculated by assuming the continuity of fs,ft,ft and 
their first derivatives at x = 0. This leads to the re- 
sults shown in Fig. 2, obtained for = z (non-colinear 
case). At zero-energy and x — > —oo, ft is dominant 
due to 7t ^ 7,/, (Fig. 2, left panel) (l8l|. I now com- 
ment the behavior of the different types of correlations 
for a; > (Fig. 2, right panel). In this area, TCSl (/f ) 
propagate independently of the two other components. 



with the characteristic vector fco. At e = 0, fco is small, 
so that TCSl propagate on a relatively long distance: 
this is why TCSl are usually called "long-range" corre- 
lations. However, this is true only at low energies. In- 
deed, for e = ±E^lf , the real part of fco is larger, so 
that TCSl decay more quickly. Conversely, for a; > 
and e = 0, the propagation of singlet correlations (/s) 
and TCSO (ff) is short-range. Indeed, fs and show 
a damped oscillatory behavior ruled by the propagation 
vectors k„, which correspond to a scale (hD / E^l^)^/'^ for 
Kl^ > F and £ = 0. In contrast, for e = ±E^f^, the 
fs and components do not oscillate and decay more 
slowly. Hence, it is the singlet correlations and TCSO 
which are long-range for e = ±E^[f . As a result, suffi- 
ciently far from the superconductor, TCSl become dom- 
inant at £ = 0, whereas singlet correlations and TCSO 
become dominant for e — ztE^^-^ . In other terms, one 
obtains v{e = 0) ~ a; and v{e — ±Ell^) ~ z. This leads 
to characteristic features in the energy dependence of the 
DOS, shown by Fig. 3a. For x sufficiently large, one ob- 
tains characteristic dips at £ ~ ±E^lf , due to singlet cor- 
relations and TCSO. In contrast, TCSl produce a zero- 
energy peak which reaches a maximum higher than the 
normal-state DOS iVo. Importantly, N{0) > Nq is not 
specific to TCSl. Indeed, an enhanced zero-energy DOS 
can also be due for instance to interaction effects [13i]. 
or to TCSO 15-11], as observed experimentally in S/F 
bilayersfll. However, in the present geometry, TCSl can 
be detected unambiguously due to the additional features 
discussed below. 

Due to the peculiar energy dependence of fco, the zero- 
energy DOS peak is surrounded by low-energy dips, and 
this ensemble shows a characteristic spatial scaling be- 
havior when X increases (see Fig. 3a). Let us note zLsd 
the position of the low-energy dips. In the limit of a 
vanishing F, Sd scales with a longitudinal Thouless en- 
ergy Exh = hD/x^, provided the finite energy dips are 
well separated from the low-energy peak, i.e. Ed -C E^l^ . 
When r is finite, a scaling behavior can persist. For in- 
stance, in the hmit £d <^ F, Eq. (j4|) gives 



TETh/2 (X l/x 



(5) 



The scaling behavior of the low energy DOS features rep- 
resents a " smoking-gun" for the fact that some supercon- 
ducting correlations propagate along the strip with the 
vector fco. Importantly, the DOS dips at £ ~ ^E^l^ fol- 
low an equivalent scaling behavior, due to the structure of 
the fco- vectors. It is instructive to discuss other magnetic 
configurations. Figure 3.b shows the DOS of the NM in 
the colinear case. The finite energy dips are still present, 
but there is no low-energy features because TCSl are ab- 
sent. At last. Fig 3c. shows the DOS of the NM when 
there is a FI contact with a uniform magnetization at 
a; < 0, but no FI contact or E^[^ = at a; > 0. In this 
case, for a; > 0, weak DOS dips appear at £ ~ ±E^l^ , 
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FIG. 4: Predictions for fhn — z, a. value of 7t much larger 
than in Figs. 2 and 3, and finite values of di^i^jiy Panel a 
shows \ft\ and \ft\ at zero energy, as a function of the 
coordinate x. Panel b shows 5N as a function of e and x. 
Panel c shows the dependence of Sd on x, for various values of 
r. The value F — 0.02A is used in panels a and b. All curves 
have been obtained with the numerical approach, except the 
dotted lines in panel c, which correspond to Eq.((5]). 

due to the spin-split minigap effect occurring for x < 0. 
Importantly, these dips are very different from those ob- 
tained in the previous cases: they are not surrounded 
by peaks and quickly vanish with increasing x. In con- 
trast, there is still a spatially-scaling zero-energy peak, 
although there is no TCSl in the circuit. This is because 
singlet correlations and TCSO propagate with the char- 
acteristic vector ko in this case. One can conclude that it 
is important to see DOS dips appear at e ~ '^^ll^ ^nd 
persist with increasing to confirm the existence of E^l^ 
at the right side of the strip. In this case, the spatially- 
scaling zero-energy peak can only be due to TCSl, which 
is the only component which can propagate with fco . For 
a good visibility of the finite-energy DOS dips, it is im- 
portant to use Ell^ < A. In practice, this limit can be 
reached by using an appropriate value for h. Note that 
the long-range behavior of TCSO and singlet correlations 
at e = iE^l^ is usually not discussed for standard fer- 
romagnets, in which exchange fields are too high. 

The amplitude of the DOS variations in Figs. 2 and 3 
is very weak due to the small 7t used. Experimental ob- 
servations call for an increase in jt- This requires to use 
a numerical approach to solve the problem in the non- 
linear limit. In practice, ferromagnetic domains have a 
finite size, so that cIl and da must be finite. The re- 
sulting finite size effects are studied below with the nu- 
merical approach. Figure 4 presents results for the non- 
colinear case. The spatial behaviors of /f , and fs 



remain qualitatively similar (see Fig 4. a). At x > 0, 
still yields the expected low-energy features in the 
DOS (Fig. 4b). The variations of Ed with x are shown 
in Fig. 4c for various values of F (symbols). The scal- 
ing behavior of the DOS low-energy features is robust to 
finite size effects. One can check that the semi- infinite 
expression Eq. (O yields the right order of magnitude 
for Ed at the right boundary of the NM strip. Finite 
size effects slow down the scaling behavior of the low- 
energy features at x — y d^, i.e. one has ded/dx — > 0. 
Nevertheless, it is still possible to observe a strong de- 
crease of Ed with X, when x is not too close to da. If 
7t is too large with respect to 7^, the zero-energy DOS 
peak has a reduced amplitude because /f is weak[l8j|. 
To maximize -/V(O), one must use 7t = 7^ and decrease 
r. For instance, using 7t = 0.08, F = 0.005A and the 
other parameters of Fig.4.b (in particular, 7^ — 0.08), 
one obtains N{0) = 1.35-/Vo at x — 40/i. Nevertheless, 
using unfavorable parameters like those of Fig.4.b, one 
still obtains N{0) = 1.00325iVo at x = 40/i, and ampli- 
tude which is measurable, in principle The geometry 
discussed here also presents the advantage that spin-flip 
scattering effects, which could reduce the amplitude of 
odd-frequency superconducting correlations, are usually 
small in NMs. Interestingly, a NM strip with a S contact 
but no ferromagnetic contacts has been studied experi- 
mentally [Ell^ = for any x). In this case, one can 
only have singlet correlations propagating with fco- As 
a result, a spatially-scaling zero-energy dip has been ob- 
served in the DOS of the NM[i,[3H. 

To conclude, TCSl can appear in a diffusive NM strip 
in contact with a S and a FI with several non-colincar 
magnetic domains. These correlations induce in the DOS 
of the NM a low-energy peak surrounded by dips, which 
show a characteristic spatial scaling behavior away from 
the S contact. Meanwhile, if the thickness of the strip 
is chosen properly, superconducting correlations between 
opposite spins will focus at finite but subgap energies. 
The simultaneous observation of both effects would en- 
able an unambiguous identification of TCSl. 

I thank T. Kontos and W. Belzig for useful discussions. 
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